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Abstract 



We propose a new mechanism for leptogenesis, which is naturally realized in models 
with a flavor symmetry based on the discrete group ^4, where the symmetry breaking 
parameter also controls the Majorana masses for the heavy right handed (RH) neutri- 
nos. During the early universe, for T > TeV, part of the symmetry is restored, due 
to finite temperature contributions, and the RH neutrinos remain massless and can 
be produced in thermal equilibrium even at temperatures well below the most con- 
servative gravitino bounds. Below this temperature the phase transition occurs and 
they become massive, decaying out of equilibrium and producing the necessary lepton 
asymmetry. Unless the symmetry is broken explicitly by Planck-suppressed terms, the 
domain walls generated by the symmetry breaking survive till the quark-hadron phase 
transition, where they disappear due to a small energy splitting between the vacua 
caused by the QCD anomaly. 
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1 Introduction 



Leptogenesis provides a natural scenario to explain the observed baryon abundance, where 
an asymmetry in the leptonic sector is initially generated by the out-of-equilibrium decay of 
heavy right-handed (RH) neutrinos and then distributed into baryons by so-called sphaleron 
processes [1]. The most attractive feature of leptogenesis is that it relies uniquely on the 
addition to the Standard Model of three RH neutrinos, addition which is already very wel- 
come to explain the smallness of light neutrino masses via the see-saw mechanism [2]. This 
synergy between neutrino physics and the origin of matter underlines the need for studying 
leptogenesis in models with realistic neutrino mixing. 

While our understanding of the origin of fermion masses and mixing angles remains at a 
primitive level, flavor symmetries provide a practical mean to replicate, within experimental 
uncertainties, the observed data. As originally proposed by Froggat and Nielsen [3], a broken 
U{1) flavor symmetry can be responsible for the small ratios between masses in the quark 
sector. Similarly, tri-bimaximal neutrino mixing, in excellent agreement with observations 
from neutrino oscillation experiments, can be reproduced within models with broken discrete 
non-abelian symmetries in the leptonic sector where sub-leading corrections in the symmetry 
breaking parameter account for small deviations from exact tri-bimaximal mixing [3]. 

In this article, we study leptogenesis in models with natural tri-bimaximal mixing from 
a flavor symmetry based on the discrete group A4 [5]. The arguments we use, however, 
can be easily extended to a more general class of models where the scale of lepton number 
violation is introduced by a spontaneous symmetry breaking mechanism. Symmetries that 
are spontaneously broken today might be restored during the early universe, due to finite 
temperature effects. We find that this symmetry restoration provides a natural scenario for 
thermal leptogenesis to work efficiently even at low reheat temperatures TfjH ^ 10^^^ GeV, 
well below the most stringent bounds from gravitino overproduction, with no need for the 
addition of structure or fine-tuning. Indeed, a crucial aspect of the y44-based models is 
the structure of the breaking sector and the particular vacuum alignment that reproduces 
tri-bimaximal mixing in the neutrino sector. The symmetry is partially broken (down to 
a subgroup isomorphic to Z2) by the vacuum expectation value (VEV) of a fiat direction 
that obtains a potential only in the presence of soft supersymmetry (SUSY) breaking terms. 
During the early universe, the finite energy density of the infiationary vacuum or the differ- 
ent occupation numbers of fermions and bosons in the thermal bath after reheating, break 
SUSY and contribute a (positiv^ mass-squared term for flat directions [9], restoring part 

^For particular forms of the Kahler potential is it also possible to arrange for a negative mass of the 
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of the flavor symmetry. The relevant aspect for leptogenesis is that, in these models, the 
lepton number violating Majorana masses of RH neutrinos are also proportional to the 
flavor symmetry breaking parameter and therefore vanish in the early universe. Thus, at 
early times, the would-be heavy RH neutrinos are effectively massless states and remain in 
equilibrium with the thermal bath down to temperatures of order the SUSY breaking soft 
masses, rh ~ 0(10^ — 10^) GeV. Below this temperature, the thermal corrections to the 
flat direction potential become too small and its origin unstable: a smooth phase transi- 
tion takes place in which the flat direction VEV rolls to large values, breaking the flavor 
group and lepton number and giving large Majorana masses Mj to the RH neutrinos, which 
suddenly flnd themselves out of equilibrium at T ^ Afj. This is the perfect starting point 
for standard thermal Leptogenesis: an equilibrium abundance of RH neutrinos, ready to 
decay via CP-violating interactions producing a large lepton asymmetry. At temperatures 
above ~ 100 GeV, sphaleron processes convert this asymmetry into the observed baryon 
asymmetry [TO] . 

A potential problem of discrete symmetries broken spontaneously at low temperatures is 
the creation of domain walls. Their energy density red-shifts slower than that of matter or 
radiation and eventually, independently of their initial abundance, they come to dominate 
the energy density of the whole universe, leading to cosmological consequences incompatible 
with observations, such as imprinting large signatures in the cosmic microwave background 
We will show that, in the /l4-based models, the putative discrete symmetry is not one 
in fact: it is broken explicitly by the QCD anomaly that lifts the degeneracy between vacua 
separated by domain wall and drives them to collapse. Thus, unless the symmetry is also 
broken explicitly by the gravitational interactions (in which case the domain walls do not 
form at all), the domain wall network disappears at the QCD scale and standard cosmology 
is recovered. 

In the next section we review models based on the A4 flavor symmetry to reproduce 
tri-bimaximal mixing, paying particular attention to the structure of the symmetry breaking 
sector. In section |3] we discuss the dynamics of the flavor symmetry breaking during the 
early universe and the thermal production of RH neutrinos. In section |4] and [5] we conclude 
with a study of leptogenesis and a note on the domain wall problem. 

flat-direction during inflation. In this case, the pecuhar dynamics of flat directions in the early universe 
[6j [7] can reproduce the observed baryon asymmetry in a non-thermal way, as in [7] or via Affleck-Dine 
baryogenesis [8]. 
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2 Tri-bimaximal Neutrino Mixing from a Discrete Sym- 
metry, at T = 



Fits to neutrino oscillation data [12] give the following values for the mixing angles, 

sin^ e,, = 0.01+°;°}^ sin^ 9,, = 0.50t',Z sin' 6,, = 0.304+n?6, (1) 

which is compatible with tri-bimaximal mixing [T3] . 

sin^ = 0, sin' 023 = 1/2, sin' On = 1/3. (2) 

It is tempting to explain this peculiar mixing pattern in terms of symmetries. In [5], a 
model based on the discrete group A4 was proposed, providing a dynamical mechanism 
which naturally reproduces this particular mixing matrix. A crucial role in this picture is 
played by the sector responsible for breaking the A4 symmetry: only for a particular vacuum 
structure can the tri-bimaximal form of the mixing matrix be reproduced. In what follows 
we will review how this works. 

The full flavor symmetry group is ^4 x [/(I) fn x Z^, where the part gives the required 
form of the mixing matrix, while the U{1)fn Froggatt-Nielsen factor provides the usual 
hierarchy among Yukawas for the charged leptons, through insertions of VEVs of the fiavon 
fields 9, 9 suppressed by some unspecified large scale A. The extra Z3 guarantees that 
charged leptons and neutrinos couple to two different set of fields at leading order. Thank 
to this separation, the symmetry can be broken down to two different subgroups in 
the neutrino and charged lepton sectors; the mismatch between the tw different residual 
symmetries is what gives origin to the mxing. The flavor symmetry breaking sector includes 
two A4 triplets, ip^ and ip'^ and two singlets, ^ and ^ which couple directly to the lepton 
sector, but also a set of "driving fields", ip^, tp^ and ^0 that are needed to build a non- 
trivial scalar potential in the symmetry breaking sector. Beside the usual field content of the 
MSSM, there are also three RH neutrino fields, in a triplet of A4. The charge assignments 
are as follow Jf| (see [5] for a review on the group) : 
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•^The U{1)fn and Z3 charge assignments given here are different from those of references [SJH]. This 
choices doesn't alter their conclusion, but will be relevant for our discussion of leptogenesis in section |4] and 
domain walls in section [Sj Note furthermore that, while the A4 factor is crucial to obtain tri-bimaximal 
mixing, the can be substituted by any (larger) symmetry that keeps the charged and neutrino sector 
separated. 
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where included is also a f/(l)_R symmetry, containing the usual R- Parity of SUSY, that 
differentiates between matter fields, symmetry breaking sector fields and driving fields. The 
superpotential can be divided in two parts, 

W = Wi + Wd. (3) 

The lepton part includes, at leading order, the terms responsible for lepton masses and is 
given by (for details on how the indices are contracted, see [5]) 

T T 

T 

+ yrT\^l)"H, + y e\NH)H^ + {xAi + iAi){N'^N-) + xb^^^^N^N') + ... (4) 

where A is the cut-off of the theory and t = (0)/A ^ (0)/A is the Froggatt-Nielsen factor 
that reproduces the mass hierarchy in the charged lepton sector. We shall absorbe the factors 
of t in a redefinition of the yukawas, 

Vet^ = 2/e, Vi^t = Vi,, = y, (5) 

which allows us to treat y as a. free parameter in the following discussion. The driving part, 
which produces the right vacuum alignment, is 

+ gi{i^l^''^^) + g2i{vl^'') + gM'^''^^) + g^ioi'' + ^75^0^^ + ^?6eof' + ... (6) 

Here the coupling constants xa,b are generally complex, while the Yukawas y,ye.n,T can always 
be made real by a redefinition of the fields; dots stand for higher order operators in the 1/A 
expansion. The mass-parameter M is taken to originate microscopically. A crucial aspect 
of this model is the existence of a fiat direction at all orders in perturbation theory (in the 
limit of exact SUSY). This can be seen as follows. Since the fields in the fiavor symmetry 
breaking and driving sectors are singlets under the gauge group, in the supersymmetric limit 
the scalar potential involves only F-terms (assuming canonical Kahler terms). 



dW ^ 



with 0i = <^Jif,^(o),^. (7) 



This means that the conditions for the minimum of the potential are simply 

dW _dW _ . _ . . . 
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and 

aw aw aw 

where the subscript i indicates the /I4- component of the triplets v^j-q'^ . Furthermore, since the 
driving fields carry two units of i?-charge, the superpotential must be linear in these fields 
(see eq.Q, at the renormalizable level). This means that the 8 equations with 7 variables 
expressed in eq.Q always have the trivial solution {^Q'f) = (^0) = 0. On the other hand, the 
remaining 7 equations in 8 variables expressed in eq (jH]) do not contain the driving fields and 
do not exhibit a unique solution, but rather a flat direction of solutions, undetermined by 
the minimum condition. This argument depends only on the particular field content and the 
charge assignments and is valid also when higher order (in A~^) terms are included. Given 
the superpotnetial Wd, eq. (EI), the conditions of eq. can be solved to find the vacuum. 
A possible choice for the VEVs of the relevant fields at the minimum reads 

ii) = 

iO = n 
(^^) = (^;^,0,0) 

(V^^) = {v',v',v') (10) 

with v'^ = —3M/2g, while u and {v^Y = —{gA/^gsju"^ remain undetermined, parametrizing 
the degree of freedom along the fiat direction (from now on we shall call u this flat directioij^. 
As explained above, in the limit of unbroken SUSY, the VEVs of the driving fields vanish 
exactly. An important feature of this model is that A4 is broken by the VEVs of two different 
fields: (f'^, which leaves invariant a subgroup Z3 C A4 and gives masses to the charged 
leptons, and ip^ , which breaks the Z3 but leaves a. Z2 G A4 invariant and gives Majorana 
masses to the neutrinos. The C ^4 acts on the A4 triplets as T = diag(l, cu, cj^), where 
u = e 3 . 

Supersimmetry breaking effects, when accounted for, contribute to the potential eq. (I7j), 
through soft masses for all scalars and A-type terms. Assuming that the main source of 
SUSY breaking preserves the flavor symmetry [T3], the main consequence of this breaking 



^Equivalently, in the basis 



e + W^gs/gii^f + <^f + <^f ) /4 



:C + tV^gs/gii^f -<P2~ 'pD) /4 ^^^^ 
X - ^^/ig^i vf) 1 2 



^W^g3.lgi{f2 - ))y 

u represents the flat direction, while u' = u" = u'" = 0. 
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will be that small VEVs for the driving fields are generated, of order fh. Furthermore a 
potential for the fiat direction is also generated, removing the degeneracy of the minima 
^ (^) = u. Its potential can be written in the generic form 

V{u) = '^mM^ + l%u' + h.c. + ..., (12) 
V / 2 ' ' 3 A ^ ' 

where the first term comes from SUSY breaking soft masses, while the second one is induced 
by non-renormalizable terms in the superpotential once the driving fields obtain VEVs of 
order m (for instance the soft term ^{1^^ ip^)^^'^ / A in the Lagrangian contributes to the cubic 
term once driving fields obtain their VEVs)lj. The coefficients a and b depend on the details 
of the model. We are interested in the case where the parameter a is negative. In this case, 
the true minimum of the potential is displaced from the origin u = and becomes very large, 
u = {a/b)A. This, together with the usual Higgs VEV Vu = {hu), generates Majorana and 
Dirac mass matrices in the neutrino sector: 

/I 0\ / A + 2B/3 -B/3 -B/3 \ 

= yvu \ 1 \ , M = i -B/3 2B/3 A-B/3 \u , (13) 

\ 1 / \ -5/3 A - B/3 2B/3 ) 

where A = 2xa and B = 2xbv^/u. Note that the Majorana mass for the RH neutrinos 
is proportional to the fiat direction VEV u and is in general very large. The heavy RH 
neutrinos can be integrated out and the mass matrix for the light neutrinos becomes m,^ = 



(m 



^^^M ^m^ = —v^y'^M ^. Its eigenvalues are 



™, = (14) 

where tan/? = Vu/vd and v = \/vl + ~ 174 GeV, while Vd = (Hd). The mixing matrix, 

U*M-'U^ = diag(Mr\ M3-1, M3-I), , (15) 
with Mi real and positive, has the tri-bimaximal form (up to a phase), 

/ 7173 \ 

= UtbUph, with Utb = -1/VQ 1/73 -1/V2 , (16) 

V -1/V6 1/V3 I/V2 J 

which is equivalent to eq. ([2]). Here Uph = diag(e*°^, e*"^, e*°^), encodes the phase information 
from the diagonalization of M. Small deviations from pure tri-bimaximal mixing can be 
accounted for by next-to-leading-order corrections in the expansion 1/A. 



■^Rcnormalization effects are expected to drive the negative mass parameter — |a|m^ to positive values at 
high energy or, equivalently, at large values of u, thus stabilizing the VEV of the flat direction at large values 
even in the absence of the term 
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3 Symmetry Restoration and Breaking 



In this section we study the dynamics of the flavor symmetry breaking sector during the early 
universe. As explained above, the flavor symmetry is broken in this model in two steps: first 
at the scale M by the VEV of ^p'^, which leaves a subgroup Z3 C unbrokeij^, then, by 
the VEV of a supersymmetric flat direction which acquires a potential only in the presence 
of SUSY breaking terms, breaking A4 completely. During inflation, the finite energy density 
of the vacuum breaks supersymmetry (the inflaton F- or D- components are necessary non- 
zero), inducing additional soft terms of order the Hubble constant Hj. In particular this 
generates a mass term 

+ jH]u^- (17) 

where the constant c/ depends on the shape of the Kahler potential and superpotential. For 
the simplest case of minimal couplings in the Kahler potential the coefficient is positive (for 
example, c/ = 3 in F-term inflation) [9]. Although it is possible to arrange for negative values 
of C] (for instance, introducing non-minimal interactions between the inflaton and the flat 
direction in th Kahler potential, as done in Affleck-Dine bariogenesis [8]), in the following we 
will assume it to be positive and of order unity. Consequently, during inflation, for Hj ^ m, 
the curvature of the potential is positive at the origin and the flat direction VEV is driven 
to zero: the Z3 C A4 symmetry is unbroken and the RH neutrinos are massless, since they 
obtain their masses only through the last two terms of the superpotential, in eq. (jlj). 

After slow-roll, the inflaton reaches the minimum of its potential and begins to oscillate 
coherently around it, while beginning to decay. The universe, whose energy density is dom- 
inated by the oscillations, behaves as matter dominated and the Hubble parameter starts 
decreasing. Although most of the inflaton condensate will decay much later (depending on 
its decay width), a small fraction decays already during the flrst oscillations and quickly 
generates a (subdominant) thermal bath with temperature 

Tosc^ koscirripiT^HH)^^^ , with ^o.c = f i^^^^l , (18) 

\57l^ 9*{Tosc) J 

where g^, counts the number of effectively massless degrees of freedom at a given temperature 
and T^H deflnes the reheat temperature where hot big bang cosmology begins. For the range 
of temperatures we are interested in, kosc ~ 0.4. Note that Tosc can generally be bigger 
than the reheat temperature. This fact, that a thermal bath is produced promptly after 

^For simplicity, we assume that the scale M is bigger than the inflationary scale and that, consequently, 
ip^ obtains a VEV before or during inflation. This assumption, however, has no influence on our conclusions. 
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inflation (rather than there being a phase in which H decreases but there is no radiation 
yet), gurantees a smooth transition between inflation and the subsequent reheating. 

Eventually, after this phase, the inflaton decays completely and the universe remains flUed 
with radiation. The constrains from gravitino over-production [15] tell us that the temper- 
ature at this stage (the reheat temperature, Tp^h) cannot be bigger than about 10^"^'' GeV. 
This, however, doesn't constrain the temperature Tqsc before inflaton decay, as the gravitinos 
produced at this stage are in fact diluted by the subsequent expansioijl|. 

So we see that ever since the end of slow-roll inflation, the universe contains a thermal 
bath of radiation with temperature Tqsc that decreases with the expansion. Since some of 
the particles in the thermal bath (the right-handed neutrinos, for instance) couple to the 
flelds ip^ and ^ that make up the flat direction, they induce flnite temperature corrections 
to its potential. For large values of T and small values of m, the potential reads 

Vt{u) = ^ {-\a\m^ + ctT^) \u\^ + (19) 

where dots stand for higher order terms in u. The positive coefficient ct receives contributions 
from all particles in the thermal bath that couple to the flat direction (driving flelds and RH 
neutrinos) [16] . For example the contribution coming from the RH neutrinos and sneutrinos 
is 

1 f.a..2 , 20^4 2 



<^T3:f^(^18xi + ^4j, (20) 

for simplicity, in what follows, we will assume it to be of order unity ct ~ C^(l)- 

As long as T ^ m, therefore, the flat direction acquires no expectation value, u = 
and the RH neutrinos remain massless. On the other hand, their Yukawa couplings are 
independent of the flat direction VEV, see eq. (j4]), and mediate fast interactions between the 
thermal bath and the RH neutrinos (and sneutrinos), bringing them in thermal equilibrium. 
Their number density reads 

nj^ = ^,T\ (21) 
It is useful to express the RH neutrino abundance in terms of its ratio to the entropy density 

n""^ 45 

since this quantity has only a mild intrinsic temperature-dependence from g^, counting the 
number of degrees of freedom in equilibrium at any given temperature. 



^The power of four in the exponential of the relation Tqsc oc i/^^* from eq. (fT8|) . ensures that at higher 
temperatures the expansion rate is much faster and that the gravitinos generated thermally before reheating 
at T ^ Tosc are effectively diluted by the expansion and never exceed those produced around reheating, for 
T ^ Trh < Tosc- 
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The universe then continues its expansion and the temperature decreases until it becomes 
comparable to the SUSY breaking scale, T ~ m. At this stage, the temperature contributions 
to the potential in eq. ( !T9l) become smaller than the (negative) soft mass term and the 
curvature of the potential changes sign making the origin unstable: the flat direction is free 
to evolve to the true fiavor-symmetry-breaking minimum of the potential. Note that for 
T > m the tunneling probability is way too suppressed [17], and the phase transition is most 
likely to take place classically (second order) as outlined so far, on a time scale set by the 
curvature of the potential, tpt ~ m~^. 

4 Leptogenesis 

In flEl [T9l [20] , leptogenesis in the framework of models with a symmetry was studied 
and found that very high reheat temperatures, T^h ~ 5 x 10^'^ GeV are required in order 
to produce the heavy RH neutrinos via scatterings in the thermal bath. This is in contrast 
with bounds coming from the overproduction of gravitinos in the supersymmetric context, 
as commented above. The symmetry restoration, however, drastically modifies this scenario. 
Its role for leptogenesis is clear: during inflation, infiaton oscillations and reheating, the VEV 
of the fiat direction responsible for the C symmetry breaking and for neutrino masses 
is kept in the origin by its interactions with other particles. At this stage the three RH 
neutrinos are massless and can be produced thermally even for very low reheat temperatures 
T > rh, independently of their would-be masses at T = 0. When the temperature decreases 
below about m (the exact temperature depending on the couplings of ^ and (f^ with the 
thermal bath), however, the phase transition takes place and the RH neutrinos become 
massive. Their masses, Mj ^ m, makes them completely out of equilibrium and eventually 
they decay via their CP-violating interactions producing an abundance of leptons over anti- 
leptons which is then transformed into the observed baryon asymmetry by sphalerons |10j . 
as it is usually the case in Leptogenesis. 

One can worry that the highly energetic decay products of the heavy RH neutrinos might 
overproduce gravitinos by thermal scattering. The decay products of the RH neutrinos, 
however are not highly energetic. In fact, since the classical phase transition that gives the 
RH neutrinos a mass takes place over a finite time tpt ~ m~^, it is easy to show that the 
RH neutrinos decay before attaining their full masses. Indeed, their lifetime is 



which is inversely proportional to their mass Mi ~ u, which increases during the phase 




(23) 
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transition {y denotes the Yukawa matrix in the basis where M is diagonal). Hence, the RH 
neutrinos decay when tat- < tpt or, equivalently when 

M, > Mt""' = 8nm/[y^y],,, (24) 

producing a thermal bath with temperature 



Using eq. f pi|) . it's easy to see that for Yukawas bigger than \y^y\ii > 10~^^ this temperature 
is safely below the most conservative bounds for gravitino overproduction [13], T < 10^ GeV. 
Such low temperature also ensures that wash-out (and hence also flavor effects, which are 
normally present at temperatures below 10^^ GeV, [21]) are negligible, since AL = 1 inverse 
decays and AL = 2 processes are out of equilibrium for T <S Mj, which is always the case. 

Thus, the baryo-to-entropy ratio is given by 

^ ~ —eY'^ (26) 
s 23 ^' ^ ' 

where the approximate numerical factor comes from the redistribution of lepton number to 
baryon number via sphaleron interactions, while the ratio of RH neutrinos number density 
to entropy at the time when RH neutrinos are still in equilibrium is given by eq. ( l22l) . 
The CP asymmetry e vanishes in the limit of exact flavor symmetry [19] and receives its first 
contributions from the terms 

|(t>^(/iV),), and |(t>^(/iV)J, (27) 

where the subscripts s, a specify in which way the flavor indices are contracted, see [5]. Again 
we can redefine the Yukawas y\^2 = ^"^1,2 to absorb the factors of t; note that yi^2 are of the 
same order as y above. The CP asymmetries for the decay of Ni, using eqs. ( IT31I271) and refs. 
[I9l I22] are given by, 

ei = tle{y2)—[-r— — sm(a2 - ai)/i2 + 77 sm(a3 - ai)/i3 



Stt \9 Re{y2y ' ' 3 

62 = Re(?/2)^^ ( p^l^^L sin(a2 - ai)/i2 + ^ sin(a3 - a2)/23 ) (2^ 



TT V 9 Re(y2)2 ' ' 3 

f I 2 
€3 = Re{y2f— -- sin(a3 - ai)/i3 - - sin(a3 - a2)/23 
ovr \ o o 
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Figure 1: The CP asymmetries ei (red, dashed), £2 (black, soUd) and £3 (blue, dot-dashed), divided 
by rf'Re{y2Y f^i' normal (left) and inverted (right) hierarchy as a function of the lightest neutrino mass, 
mi = mi(m3) for normal (inverted) hierarchy. We used Re{y2)^ = Re(jji)^, while Am^^j = 7712 — ml = 
7.67 X 10-5 eV^ and Am^^^ = |to§ - | = 2.46(2.37) x 10"^ eV^ for NH (IH), see [Ml- 



The flavor symmetry breaking parametecl f] = \{(f'^)\/A ^ \{if^)\/A ^ (0/-^ realistic 
models is approximately [5] 

0.0022 <7]< 0.05, (29) 

where the upper bound comes from the requirement that subleading corrections to the lepton 

mixing matrix not be too large, while the lower bound stems from the requirement that the 

tau Yukawa be small enough to justify a perturbative expansion. The function fij depends 

on the ratios between RH neutrino masses as 

M, /2 - Mf/M^ ^ ^ /I + Mf/MfW 

f = lL^ _ (1 + M^/Mf) In ' , (30) 

and in the limit of strongly hierarchical RH neutrinos Mj ^ Mj goes as fij ~ (3/2)Mj/Mj. 
In the models studied here, the masses Mj and angles are completely determined in 
terms of the solar and atmospheric mass parameters and the (unknown) mass of the lightest 
neutrino, m/. Figure 1 shows that the largest^ CP asymmetry is €2 and that, 

62 ^ f 10-2 for NH 



Re{y2)W ~ I 2 for IH ' ^^^^ 

Here we have taken Re(?/2) ~ and NH (IH) stand for normal and inverted hierarchy 

respectively. Comparing our results eqs. fl26|3ip with the observed baryon asymmetry [25] 

yWMAP ^ ^Q_g7 ^ Q_Q3) ^ (32) 



^Note that, while tp'^ being a flat direction naturally obtains a large VEV even at low scales m, the VEV 
of ip^ must be generated at a much higher scale M ^ ifi in order to produce enough CP violation (and the 
right size of the r Yukawa). The fact that |(<(5^)| ~ Kv'"^)! is not surprising , as one might assume that the 
same physics that gives rise to the non-renormalizable terms in eq. (jl2p , is responsible for the mass-term M 
as happens, for example, in composite Higgs models [53], where naturally the Higgs vev is w < A. 

^The three RH neutrinos obtain masses simultaneously and decay, almost at the same time. Therefore 
the lepton abundance will be produced mostly by the RH neutrino with the biggest e. 
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we find the following constraints on the Yukawa couplings, 

10-' <y^^ yl ^ yl < 1, (33) 

for normal hierarchy, while 

10"' <y^-yl- yl < io~', (34) 

for inverted hierarchy. The uncertainties depend on the value of f], eq. ( 129|) . In any case, 
this is largely compatible with the bounds given above for gravitino overproduction. 



5 Domain Walls 



Spontaneously broken discrete symmetries are commonly thought to be incompatible with 
standard cosmology, due to the existence of domain walls [11]. Domain walls are stable 
field configurations which interpolate between regions with different, but degenerate, vac- 
uum states. As discussed above, spontaneously broken symmetries are generally restored at 
high temperatures and the phase transition occurs only when, during the expansion of the 
universe, the temperature is lowered below the critical value. The problem is that different 
regions of size ~ are causally disconnected and, once the symmetry is broken, will find 
themselves in different vacua, divided by domain walls (the correlation length is reduced to 
in the case of a second order phase transition). The energy density of domain walls is 

a 

Pdw^ (35) 

where R is the typical size of the domain walls, while a is the surface energy density. R is 
proportional to the scale factor r and the energy density of domain walls scales as pdw ~ 
much slower than radiation prad ~ or matter, p^at ~ Therefore, if the phase 

transition takes place after inflation, these domain walls will eventually come to dominate 
the energy density of the universe leading to unacceptable observable signatures in the cosmic 
microwave background. On the other hand, if the scale of the phase transition is higher than 
the inflationary scale, then one can assume that the universe expanded enough that the 
relevant portion around the present horizon is clear from domain walls. This is the case for 
the if'^ phase transition, assuming that the scale M in eq. ([6]) is larger than the inflationary 
scale. Unfortunately, ip^ obtains a VEV at T ~ m well after inflation, and the domain 
walls produced at this stage, separating vacua with different Z3 C charge, can lead to 
cosmological catastrophe. In this section we discuss three different solutions to this problem. 

1. The discussion above assumes that the symmetry ^4 is exact. But it's not. When 
the model is extended to include quarks [26], these transform under a chiral representation 
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which is anomalous: the A4 symmetry is only an apparent symmetry of the classical action, 
broken at the quantum level. The same mechanism that gives the axion a mass, here lifts 
the degeneracy between vacua and leads to the disappearance of the domain walls 



Referring to [26] for the details of how the ^4 symmetry can be extended to quarks, we 
recall here how quarks transform under the Z3 C A4, 

Ur Ur 

dn w"^ dn (36) 

with w = e*^'^/^, while the other quarks are singlets. This transformation is chiral and it has 
a QCD anomaly, in the sense of ref. [27] - it acts non-trivially on the quark mass matrix Mq. 
Since a phase shift in argdet(Mg) corresponds to a shift in the ^QCD-angle, we see that the 
action of Z3 C A4 interpolates between 6qcd = 0, 27r/3, 47r/3. Now the same physics that 
gives the axion a mass, lifts the original degeneracy by a small amount Ap, proportional to 
the axion mass, 

Ap ^ f^ml^^ (37) 

at T = 0, where m^, f^^ are the pion mass and decay constant and rrtu^ are the m, d- 
quark masses. At large temperatures, this small energy splitting has no influence on the 
dynamics of domain walls, but at temperatures below the chiral phase transition, T < 
Tqcd = 250 MeV, the pressure caused by this energy difference will become sizable and the 
domain walls collapse, [27] . We need to make sure that this collapse takes place before the 
domain walls dominate the energy density of the universe, which would otherwise lead to 
power-law inflation that empties the universe from matter and radiation. 

In our model the domain walls have width S ~ and false vacuum energy p^ ~ iv?v?, 
resulting [TT] in a surface energy density a = 5p^ ~ mn^. The tension, with a force per unit 
area Ft ^ a/R, tends to drive small walls to collapse, while wiping away inhomogeneities on 
larger walls. On the other hand, friction Ff ~ —T^v {v being the wall velocity with respect 
to the bath with temperature T), caused by particles changing their masses while crossing 
the wall, retards this collapse and only scales much smaller than the horizon are wiped away. 
The size of the typical scales surviving is set by a balance between tension and friction, 

^ = T%. (38) 

Indeed, the time required to straighten a wall of size R is 

R R^T^ , , 

td = -^ , 39 

V a 
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which means that in one Hubble time, td = H ^ ~ MpJ {gl^'^T'^), only scales bigger than 

R - ^ (40) 

will remain. This sets the typical size of domain walls at temperature T. Using eq. fl35|) we 
can compare the energy stored in the domain walls at the QCD phase transition with the 
energy density of the universe at temperature Tqqu'. 



Prad ngt^^\J MpiTgcD 1010 GeV ^ ' 

In order for the domain walls never to dominate the energy density of the universe, we need 
that the flat direction VEV, 

u < 10^° GeV. (42) 

Recall from eq. (fT3l) that the fiat direction VEV is of order the RH neutrino masses. Mi ~ u, 
and that these are fixed through eq. ( fT4l) by the size of the Yukawas y and the light neutrino 
masses. Thus, the bound of eq. (142|) can be thought in our model as a bound on the Yukawa 
couplings, 

y^y < lo-^ (43) 

where we have used m; = rriatm = 0.05 eV. Comparing this with the results of eq. fl5^ . 
shows that this solution of the domain wall problem is compatible with leptogenesis only 
in the IH case, while the NH case is disfavored since there domain walls would come to 
dominate the energy density of the universe before their collapse. The fields (f^ and ^ along 
the fiat direction also carry a non vanishing Z3 charge which is liable to yield to its own 
domain walls. This symmetry however is violated explicitly by the /i-term of SUSY, fiHuHd- 
Due to this relatively large explicit violation, similarly to the mechanism just showed, the 
domain walls will collapse long before the QCD phase transitior0. 

2. It must be noted that the mechanism explained above relies on the fact that the discrete 
symmetry A4 is not violated by the gravitational interactions, which is true, for instance 
if the symmetry originates within a continuous gauge group. However, if a small, Planck- 
suppressed symmetry breaking term is present in the effective Lagrangian, the situation is 
different. Consider for example the term 

1 



M, 



pi 



(y.V')'(¥'V')eo, (44) 



^"Alternatively, the symmetry - its only role being the separation between charged and neutral lepton 
sectors - can be substituted by a larger discrete symmetry (such as Z4), with anomalous charge assignments. 
In this case its domain walls will disappear together with the ones from the broken Z3 C A4 symmetry. 
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which exphcitly breaks Z3 C A4, since the contraction (...)' is not invariant. The inclusions of 
these terms does not compromise the vacuum ahgnment, since they are suppressed by at least 
a factor A/Mpi with respect to symmetry preserving next-to-leading-order operators, which 
already give contributions of order 77, [5]. When (f^ and ^0 acquire VEVs, the contribution 
from this term would be enough to erase completely the barrier separating different vacua. 
In this case, independently on the value of u, the domain walls will not even form. Indeed, 
the energy difference between vacua provided by the term eq. (14^ is of order mu'^/Mpi, and 
for u > 10^° GeV is bigger then the height of the barrier between vacua, which is of order 

3. Before concluding, we note that in more conventional models of leptogenesis (or where 
the barion asymmetry is produced through another mechanism, such as |^ or Affleck-Dine 
bariogenesis [8]), the domain wall problem can be solved through a non-restoration of the 
discrete symmetry at high temperatures. Indeed, it is possible to arrange for (non-minimal) 
couplings in the Kahler potential between the inflaton and the fields making up the fiat 
direction (^ and (f^), such that the term of eq. (fT7|) setting the curvature of the fiat direction 
during infiation, is negative [9]. In this case a phase transition takes place during infiation 
and the domain walls produced are effectively diluted by the subsequent expansion. Normally 
the symmetry would be restored during the reheating phase, reintroducing the domain wall 
problem. However, since fiat directions acquire very large VEVs in comparison to their 
curvature, they induce very large masses to all fields which couple to them in a renormalizable 
way. These heavy fields decouple from the effective theory and are no longer able to mediate 
temperature effects to the flat direction, so that its potential is approximately given by the 
T = potential of eq. ( TT2l) . In this case the symmetry is never restored at high temperatures 
and domain walls do not cause any problem. 

To summarize, if the condition eq. fj42|) is met, then the domain walls collapse at the 
QCD phase transition without consequences. Otherwise, a small gravity mediated explicit 
breaking of the symmetry is needed. Alternatively, if the baryon asymmetry is produced 
by a different mechanism than the one discussed in the flrst part of this article, it might be 
possible to arrange for the symmetry to be broken already during inflation: in this case the 
symmetry is never restored and domain walls do not form. 

6 Conclusions 

We have analyzed thermal leptogenesis in the framework of models with an A4 flavor sym- 
metry that naturally reproduce tri-bimaximal neutrino mixing, although our flndings ap- 
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ply to the more general case where the RH neutrino masses arise from a mechanism of 
spontaneous symmetry breaking. We found that leptogenesis can be successful, without 
any fine-tuning, even at low energies well below the gravitino bound, independently of the 
masses of RH neutrinos. This is possible thank to a mechanism of symmetry restoration, 
aheady present in the y44-based models, that ensures that the RH neutrino masses vanish 
at temperatures T > 1 TeV. Hence, during the early universe, the RH neutrinos can be 
produced copiously by thermal scattering, without requiring very high temperatures that 
would lead to the overproduction of gravitinos, in contrast with the standard cosmological 
scenario. The symmetry breaking discussed here takes place along a supersymmetric flat di- 
rection, lifted only by SUSY breaking effects. Its potential is therefore naturally almost flat 
and induces very large field VEVs, despite its scale being of order the soft SUSY breaking 
masses, rh ~ O{10^ — 10^) GeV. Hence the phase transition that gives a mass to the RH 
neutrinos takes place only at temperatures much smaller than the would-be RH neutrino 
masses. Eventually the heavy neutrinos decay out of equilibrium via their CP-violating in- 
teractions producing an abundance of leptons over anti-leptons which is then transformed 
into the observed baryon asymmetry by sphaleron interactions. 

The /l4-based models are very predictive in terms of the neutrino mixing angles. Indeed 
all 3 mixing angles are predicted to depart from their pure tri-bimaximal values, eq. ([2]) 
by the same amount. Since measured deviation of 612 are very small, see eq. ([1]), the same 
must be true for 613. So, observation of ^13 7^ in future experiments will be able to exclude 
models with tri-bimaximal mixing in the neutrino sector. 

Models with discrete symmetry spontaneously broken at energies below the inflationary 
scale, generally suffer from a domain wall problem. We showed that the yl4-based models, 
when extended to account for quark masses, posses a QCD anomaly, in the sense that the 
A4 acts non-trivially on the quark mass matrix. In this case, the same mechanism that 
gives a mass to the axion, here lifts the degeneracy between vacua on opposite sides of 
the domain walls. We found that for neutrinos with inverted hierarchy this solution of the 
domain wall problem is compatible with our scenario of leptogenesis. For neutrinos with a 
normal hierarchy, on the contrary, this effect is not enough to avoid a domain wall dominated 
universe if we insist on the mechanism of leptogenesis as described above. In both cases, if 
the discrete symmetries are broken explicitly by the gravitational interactions, there remains 
no barrier separating different vacua and domain walls never form. 
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